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Abstract 

The modified perturbative approach in which transverse degrees of freedom 
as weh as Sudakov suppressions are taken into account, is apphed to B decays 
into two vr mesons. The influence of various model parameters (CKM matrix 
elements, B decay constant, mesonic wave functions) on the results as well 
as short distance corrections to the weak Hamiltonian are discussed in some 
detail. The perturbative contributions to the B decays yield branching ratios 
of the order of 10~^ — 10~^ which values are well below the upper limit for 
the B^ — > vr^vr^ branching ratio as measured by CLEO. 
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I. INTRODUCTION 



Recently the rare decay modes of the B meson into hght mesons, e.g., —>■ vr+vr", 
attracted much attention of theoreticians although such decays have barely been observed 
(see, however, ref. The reason for this interest lies in the expectation that, owing to 
the large momentum release, such processes are amenable to a perturbative treatment. If 
that is the case the rates for such B decays can be calculated in a similar fashion as for 
instance the large momentum transfer behavior of the pion's electromagnetic form factor. 
The perturbative approach would form an attractive alternative to phenomenological models 
based on soft physics (e.g., overlaps of wave functions and/or vector meson dominance 
A reliable estimate of the rates F for exclusive non-leptonic B decays is also of enormous 
importance in the investigation of CP violation. For the decays of charged B mesons, for 
instance, one studies asymmetries between the rates of CP conjugated processes. A non- 
zero rate asymmetry, signaling CP violation, requires phase differences in the interfering 
amplitudes which may be generated by penguin diagrams. 

The first perturbative calculation of exclusive B decays has been carried out by Szczepa- 
niak, Henley and Brodsky . Other authors have repeated that analysis using similar 
methods. The results for the decay widths obtained by the various authors agree with each 
other to a certain extent. However, there is a technical difficulty which has not been solved 
in a satisfactory way by these authors. As a consequence of the collinear approximation used 
in these analyses, i.e. of the neglect of the transverse degrees of freedom, a singularity ap- 
pears in the end-point regions where one of the partons carries most of the momentum of its 
parent hadron. It is presumed that a Sudakov factor damps the end-point regions sufficiently 
and a finite, stable result remains. This conjecture is taken as a justification for cutting off 
the end-point regions, a procedure which leads to an unwelcome strong dependence of the 
results on the position of the cut-off. 

Botts and Sterman have recently calculated a Sudakov factor, comprising gluonic 
radiative corrections, for exclusive reactions in next-to-leading log approximation using re- 
summation and renormalization group techniques. This approach necessitates the explicit 
consideration of the transverse degrees of freedom which, as we mentioned above, are con- 
ventionally neglected. The Botts-Sterman approach which may be termed the modified 
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perturbative approach, has aheady been apphed to several exclusive reactions, among them 
is the pion's and the nucleon's electromagnetic form factor |IU|-[T^. It has turned out that 
the modified perturbative approach allows a reliable, theoretically self-consistent calculation 
of the perturbative contributions to electromagnetic form factors. Previous objections [15 



against the applicability of the standard perturbative approach do not apply to the modified 
approach. 

The aim of the present paper is to perform a detailed analysis of B decays into two vr 
mesons using the modified perturbative approach. As it will turn out, the above-mentioned 
singularity disappears, hence there is no need for cutting off the range of integration. A 
second advantage of the modified approach is that it allows both a simple parameter-free 
treatment of the strong coupling constant and a choice of an appropriate renormalization 
scale which keeps under control contributions from higher order perturbation theory. Owing 
to these two advantages the modified perturbative approach provides a reliable estimate of 
the perturbative contributions to B decays into two tt mesons. This will be elucidated in 
some detail in the body of the text. Comparison with future experimental results will reveal 
whether the processes under investigation are indeed dominated by perturbative contribu- 
tions or whether they are still under control of soft physics. It should be emphasized that of 
all experimentally accessible exclusive decays into hadronic final states the process B vrvr 
has the largest momentum release and can therefore be considered as a hard processes best. 
If here the perturbative approach fails in comparison with experiment its application to 
other exclusive non-leptonic decay processes seems dubious. 

In Sec. 2 we discuss the calculation of the decay rate for the process B^ tt+tt^ within 
the modified perturbative approach. The numerical results and a discussion of a number of 
approximations made in the calculation are presented in Sec. 3. Short distance corrections 
will be discussed in Sec. 4 as well as a calculation of other B —>■ tttt decay modes. The weak 
B ^ n transition form factors can be calculated in a similar fashion as the B ^ tttt decay 
rates. Results for transition form factors are presented in Sect. 5. The paper terminates 
with a few concluding remarks (Sec. 6). In an appendix some details about the Sudakov 
factor are presented. 
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II. THE MODIFIED PERTURBATIVE APPROACH TO B DECAYS 



We start with the process, vr+vr , and consider the unperturbed weak Hamiltonian 

G 

Hw = 4.-^Vu{uai^,Lha){dp^^Lup) (2.1) 

where a and (3 are color labels. L is the left-handed projection operator. G^? is the usual 
Fermi constant and Vu = KTdKtfe where the Vij denote the CKM matrix elements. Because 
of the large momentum release in the decay of the B meson into two light mesons it seems 
possible to apply perturbative QCD. For the same reason one may also neglect long-range 
final state interactions. On these premises the decay amplitude M. = {71^11^ \T-lw\B^) can 
be calculated along the same lines as for instance the pion's electromagnetic form factor. 
With the help of the factorization formula of perturbative QCD for exclusive reactions the 
decay amplitude is written as a convolution of meson wave functions and a hard scattering 
amplitude Th to be calculated from the effective Hamiltonian (|2.1|) and, to leading order 
in the strong coupling constant as, the exchange of a gluon between the spectator quark 
and one of the other quarks (see Figs. 1 and 2). In principle there are contributions from 
higher Fock states too where additional quarks and/or gluons appear in the mesons. Such 
contributions will be ignored by us since they are suppressed by powers of as/M"^ {Mb 
denotes the mass of the B meson). In order to calculate the decay amplitude we utilize 
the modified perturbative approach. The crucial element for incorporating the Sudakov 
corrections in the calculation of that amplitude within the modified perturbative approach 
is the explicit dependence on the transverse degrees of freedoms in the convolution of wave 
functions and hard scattering amplitude. This convolution formula can formally be derived 
by using the methods described in detail by Botts and Sterman 0. Adapting Li and Yu's 
result for the semi-leptonic B —>■ n transition matrix element [|l^] to our case of the B —>■ nir 
decay amplitude, the convolution formula can be written in the form 

G, 

=u.. I Id.rl 

477 

xexp[-Si{x},{b},MB)]. (2.2) 



G f 
M = J [dx] 



"^'^ v|/;^(xi, bi) ^1^{X2, bs) Th{{x}, {b}, Mb)'^b{x, -b) 



X {xi, i = 1,2) denotes the usual fraction of the Pb {Pi, P2 ) component of the B{ji = 
7r+, 7[2 = vr^) meson's momentum the quark carries. The antiquarks carry the fractions 



1 — X and 1 — Xi respectively, b (bj) is the quark- antiquark transverse separation in the 
B (vTj) meson, [dx] is short for the product dxdxidx2 and {x} stands for the set of variables 
X, Xi,X2; analogous definitions are used for the b variables. Last = B, vTj) is the Fourier 
transform of the momentum space wave function "^h- 

^^hixhM) = (^) V d^k±h^h{xh,\^±h)e-'''^'^-'"^ (2.3) 

where the Fourier transform variable k^/j is the quark's transverse momentum defined with 
respect to the momentum of the meson h. Strictly speaking the meson wave function 
represents only the soft part, i. e. the full wave function with the perturbative tail removed 
from it. The full meson state is described by the product of the (scalar) wave function 
a color and a flavor function of obvious form and of the spin wave function which is written 
in a covariant fashion 

r. = ^ih + M,) [l + ^ ]^ 75. (2.4) 

Kf^ is the relative momentum of quark and antiquark forming the meson h and is the 
meson's mass. Note that the relative momentum is only determined up to a multiple of the 
hadron's momentum p^. The term ~ ^ in ( |2.4| ) takes into account the fact that quark and 
antiquark do not move coUinear with their parent hadron; it represents the first term of an 
expansion over powers of K. The function gh{xh) is controlled by soft physics as is the wave 
function "^h, it cannot be calculated from first principles at the time being. Therefore, one 
has to rely on models if its explicit form is needed. Fortunately, as will become clear in a 
moment, we have not to specify to the order we are working. 

In ( |2.2| ) Th represents the Fourier transform of the hard scattering amplitude. To lowest 
order perturbation theory, it is to be calculated from the Feynman graphs shown in Fig. |^. 
In order to perform that calculation it is convenient to work in the rest frame of the B meson 
where the hadronic momenta in light-cone coordinates read 

P^ = (Mb/v^,Mb/v^,Ox), K = (Mb/V2,0,0^), p^=p^-p^ (2.5) 

The quark momenta are specified in Fig. |I|. In a possibly crude approximation, yet well in 
line with the basic ideas of the parton model, we chose the relative momenta Kj^ in such a 
way that they have only transverse components k_|_/i and, hence, the Kf^ are all orthogonal to 
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the mesonic momenta. The wave function for the B meson is known to have a pronounced 
peak at X = Xq = rrih/MB {mf, being the 6-quark mass). It is clear therefore that only 
regions close to the peak position contribute to any degree of significance. This implies, to 
a very good approximation, that the 6-quark mass equals xMb- In accord with the heavy 
quark effective theory it also implies equal velocities of the 6-quark and of the B meson up 
to corrections of order k^/Ms- Since the r.m.s. transverse momentum is of the order of 
300 to 400 MeV these corrections can safely be ignored. The masses of the light quarks and 
mesons are neglected in the kinematics. 

Given these assumptions and approximations the denominators of the internal quark and 
gluon propagators in the graphs Fig. 2a and 2b read 

Db = ql -ml = {l-2x + Xi) M| - 
D, = ql = -{l-x)Ml-\<il 

ql = -(1 - x)(l - xi) Ml - (kx - k^i)2 (2.6) 



where, following previous authors, e.g., we neglect terms ~ (1— x)^. If one would ignore 
the (1 — x)-term in Di, which is equivalent to the assumption = Mb-, the denominators 
of the internal partons are of the same type as in the case of the pion's electromagnetic 
form factor with the i?-meson mass playing the role of the momentum transfer. It has 
been shown 1 13,12] that for momentum transfers of the order of the value of the 5-meson 



mass perturbation theory can be applied self-consistently. In our numerical evaluations of 
the decay width we, however, take into account the (1 — x)-term in D^. It leads to a pole 
within the range of integration which is handled in the usual way by using the prescription 
l/{x — ie) = P(l/x) + m5{x). This pole which corresponds to the situation of the b-quark 
propagator going on-shell, is not a pinched singularity and is, therefore, not associated with 
long distance propagations of the b quark. Hence, a perturbative treatment of the pole 
contribution is justified |T^. For a detailed discussion of this pole we refer to p8|] . 

The hard scattering amplitude Th can now be easily worked out as traces of products 
of the mesonic spin wave functions and spinor expressions representing the graphs shown in 
Fig. 1^. From the color structure of the operator ( p.l| ) it is clear that the graphs shown in 
Figs. 2c and 2d do not contribute to the process B^ — * tt+tt^. For graph 2a, one finds 
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n ~ Tr {r^,7^L(^b+ nibh^TBY} Tr {r^,7^^} (2-7) 
where, as usual, f = 7or'''7o. Working out the traces we find the simple results 



T,({x}, {k^}) = 16v^7ra,(/i)^^^-^4^ [1 + 0{kjMs. k^i/Ms)] (2.8) 
and, since we neglect the pion mass, for the graph 2b 

1^ = 0{k^/MB, k^i/Ms). (2.9) 

\-Lb\ 

fi, the renormalization scale, will be chosen subsequently. A color factor is not included in 
Tft. The Fourier transform of reads 

4^2 



Tt{{x}, {b}) = -^as{fi)MU2x - xi)Ko{J{l - x)(l - a;i)Ms6) 



xKo{y/2x - xi - IMfilb + bi|)5(b2) (2.10) 

where ii^o is the modified Bessel function of order zero. Note that for 1 + xi > 2x the 
second Bessel function in ( p.lO|) has an imaginary argument. In this case we use the relation 



Ko{iz) = -i/2H^^\z) where H^^'^ is a Hankel function. This treatment of the Fourier 
transform of Db is in accord with the ie prescription mentioned earlier. 



The last item in j^^ ) to be specified, is the Sudakov factor exp[— S*] which comprises 
those gluonic radiative corrections not taken into account by the usual QCD evolution. On 
the basis of previous work by Collins et al. ||TU[, Botts and Sterman |^ have calculated the 
Sudakov factor using resummation techniques and having recourse to the renormalization 
group. Adapting their result to the case at hand, B rnr, we write 

S'(a;i,X2, 61,62, /i) = S^^{xi,bi, Mb, f^) + S'^2(x2, 62, M^, /i) (2.11) 

S^X^,, Mb, /i) = s(x„ 6„ Mb) + s{l - x„ 6„ Mb) - A/Po In f i^TTTT^^l (2-12) 

\\n{l/{biAQCD))J 

where jSo = 11 — 2r7,//3; Uf is the number of the active fiavors, which is 4 in our case. The 
last term arises from the application of the renormalization group. The Sudakov function s, 
which includes leading and next-to-leading logarithms, is given explicitly in the appendix. 



Contrary to |T6| we do not allow for a Sudakov factor for the B meson. We present the argu- 



ments for the ansatz ( p.ll| ) in the appendix. Although our treatment of the B meson differs 



from that of [|T6|, the numerical results are practically the same since the Sudakov factor 



for the B meson used in |T^, only yields a tiny additional suppression of the perturbative 
contribution of the order of a few percent. 



Following other authors |]10| , p!6| we choose as the renormalization scale the largest mass 
scale appearing in the process 



^/(l -x)(1-xi)Mb, 1/6, 1/61 j. (2.13) 



fx = max 

The evaluation of the —>■ tt+tt^ amplitude requires the integration over the transverse 
separations. This implies two angle integrations, say the integration over the directions of 
b, simply providing a factor of 271, and the integration over the relative angle (f between 
b and bi. The latter integration is non-trivial since the relative angle also appears in the 
argument of one of the Bessel functions (see ( p.lO|) ) but it can be carried out analytically 
by means of Graf's theorem 

f{y,bM)^^JdipKo{y\hTh,\) 

= 9(6 - h)Ko{yb)Io{yh) + Q{h - h)K^{yhi)h{yh) (2.14) 

where 6 is the step function. /o(x) = Jo{ix), where Jo is the Bessel function of order zero. 
As an inspection of ( |2.1CI| ) and ( p.ll| ) reveals the X2 and b2 integrals factorize. With the help 



of the 5-function appearing in (|2.1CI|) the b2 integration can directly be carried out, yielding 
a factor 



1 
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dx2K(x2,0). (2.15) 



The value of this integral is fixed by the vr^ — > /i^z/^ decay and equals /^/ (2v^) |2^ where 
(= 130.7 MeV) is the pion decay constant. The fact that the 7r2-part in ( p.2|) factorizes 
implies the factorization of non-leptonic decay amplitudes into a product of two current 
matrix elements 

M = ^^;.(7^-(p2)|J^M0)(7r+(pl)|J^I^°M) (2-16) 
with the first matrix element being parameterized in the usual way as 

(7^+(p2)|J(^.|0) = ^p^. (2.17) 



Frequently (p.l6|) is used as a hypothesis in soft models for decay amplitudes. 



Putting all together we have the following representation of the tt+tt amplitude 

M = -^^^^v^M%n, (2.18) 

where 

f-K f 

^b= 2^^^ J da;dxiW66id6i^;^(a;i,bi)^B(a;, -b)exp [-S^^{xi,bi, Mb, fi)] 

xa,(/i)[2x - xi]iro(V(l - x)il - Xi)Mb6)/(V2x - Xi - 1Mb, b, h) (2.19) 

The color factor is \/?>Cp where Cp{= 4/3) is the Casimir operator of the fundamental 
representation of S'f/(3)c. The remaining four dimensional integral has to be carried out 
numerically. Before presenting the numerical results (see Sec. 3) two remarks are in order: 

i) Owing to the behavior of the Kq function the hard scattering amplitude has a singularity 
of the type ln(l — a;i) for xi ^ 1, x fixed. Since the pion wave function provides a factor 
1 — a;i in that limit (see Sec. 3) the integral in ( ^.19| ) is regular. This result appears as a 
consequence of the fact that the transverse degrees of freedom are explicitly considered in 
the modified perturbative approach. In the standard approach on the other hand, the 
neglect of the transverse momenta (and the assumption mf, = Mb) leads to a 1/(1 — XiY 
behaviour of the hard scattering amplitude (see (|2.6|) ). Hence the decay amplitude is log- 
arithmically singular. That property of M. forced previous authors to cut off the xi-integral 
at Xi = 1 — e (with e being chosen to lie in the range 0.05-0.1) under the argument that the 
Sudakov factor will suppress the end-point region. It should be clear from our discussion 
that the singularity of the amplitude is avoided when the transverse degrees of freedom are 
taken into account. The Sudakov factor provides only further suppression of the end-point 
regions. In contrast to our results are therefore insensitive to the end-point regions. 
Moreover, they do not depend on a cut-off parameter. 

ii) For b{bi)AQCD —>■ 1, «s(/u) is singular. However, as an inspection of ( [Al| ) reveals, the 



Sudakov factor tends to zero faster than any power of ln{l/b{bi)AQCD) compensating the 



singularity and rendering the integral in ( 2.19 ) finite. This is an attractive feature the mod 



ified perturbative approach possesses; it allows to choose the renormalization scale in such a 
way that large logs from higher order perturbation theory are avoided and a finite result for 
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the amplitude A4 is nevertheless obtained without introducing an external parameter (such 
as e.g. a gluon mass) to regularize a^. We will take up this point again. 



III. NUMERICAL EVALUATION OF THE vr+vr DECAY WIDTH 

The first step in the evaluation of that decay width is to specify the wave functions. 
The pion wave function is rather well-known since it is constrained by the decay processes 
TT^ ^^Vfji and 7r° 77 pO|. In [jl4[ a parameterization of that wave function is given 



which respects these two constraints. The soft transverse configuration space wave function 
is written as 

il^,{xi, b,) = ^ ^^.(xi) %X^Xi{l-Xi) b,). (3.1) 

For the distribution amplitude $ we use 

$;J^(x,) =6x,(l-x,) (3.2) 

and, as an alternative, we use a distribution amplitude originally proposed by Chernyak and 
Zhitnitsky |^ 



<l>^f (a;,) = 30a;,(l - Xi) (2x, - 1)^. (3.3) 
The 6-dependent part S is assumed to be a simple Gaussian 



S^^(Jxi(l - Xi) hi) = An exp (-a;i(l - Xi) ki/AaJ . (3.4) 



More complicated forms than (|3.4| ) (e. g., a two- humped shape of the momentum space 
wave function) are proposed in on the basis of dispersion relations and duality. At large 
transverse momentum, however, the soft momentum space wave function should behave like 



a Gaussian p2[. The examination of a number of examples corroborates our expectation 
that forms of S^ri other than ( p.4| ) will not change the results and the conclusions presented in 
our paper markedly. The parameter a^r appearing in the Gaussian has a value of 2.02 GeV^^ 
for both the distribution amplitudes. 



In [0 it is also shown that the asymptotic (AS) distribution amplitude combined with 



the Gaussian (|3.4|) leads to a reasonable description of the vr — > 7 transition form factor 
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within the modified perturbative approach while the CZ distribution amphtude seems to be 
in conflict with the data. We nevertheless retain that distribution amplitude since we cannot 
exclude it deflnitely and in order to examine the influence of the form of the distribution 
amplitude on the flnal results. 

For the B meson we take the Bauer-Stech-Wirbel wave function which has been 
proven to be useful in weak decays. Writing that wave function in a fashion similar to ( p.lD , 
we have 



Es = Att exp [-byAa%) (3.5) 

and 

^b{x) = Ax{1- x) exp[-a|M|(a; - xof]. (3.6) 

The distribution amplitude $b exhibits a pronounced peak, its position is approximately 
at Xo = rrih/MB = 0.93 for nib = 4.9 GeV. This property of the 5-meson wave function 
parallels the theoretical expected and experimentally conflrmed behaviour of heavy meson 
fragmentation functions. The two parameters appearing in the 5-meson wave function, 
namely A and as, are flxed by normalizing the wave function to unity - the neglect of higher 
Fock states seems to be a reasonable assumption for heavy hadrons - and by taking a value 
of 180 MeV for the 5- meson decay constant Jb- This value has been found in a recent 
lattice gauge theory analyses Under these conditions A and qb have the values 63.05 



and 1.491 GeV~^ respectively. One may be tempted to use a non-factorizing form for 
similarly to that of in ( p.4|) . However, we favor the function ( p.5|) since the non-factorizing 
form has theoretical deflciencies in the formal limit Mb oo and is, therefore, in conflict 
with the heavy quark effective theory as has been shown in [p^ . 



Using a value of 0.975 for the CKM matrix element Vud and 200 MeV for Aqcd, we flnd 
the following numerical results for the decay width of the process n^n^ when either 

the AS distribution amlitude (|3.2| ) or the CZ one (|3.3|) is employed: 

r(S° ^ vr+vr-) = -^\M\' = ( X IQ-^V | ^'^^ ^^^^ (3.7) 
^ ^ 16vrM^' ' U.005; | 2.58 (CZ) 

Since the CKM matrix element Vub is still poorly known we tentatively take for it the value 



0.005 and show Vub explicitly in (|3.7|) . We repeat that we favor the result obtained with the 
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asymptotic distribution amplitude and believe that the CZ distribution amplitude perhaps 
leads to an overestimate of the perturbative contribution to the tt+tt^ decay. The 

predictions for the decay width are also subject to errors due to uncertainties in the model 
parameters and assumptions leaving aside the difference between the results for the two 
possibilities for the pion's distribution amplitude. The influence of small modifications in 
the wave functions like the use of the BHL factor exp[— a|mV(a;j(l — Sj))] in the distribution 



amplitudes (where ruq is a constituent mass) |2y] or changes in the b dependence is rather 
mild, the decay width is typically altered by 5 — 10%. The largest uncertainty in our results 
arises from the error of 50 MeV for quoted in [^]. Although this error is very large and, 
moreover, the decay width depends quadratically on fs , the resulting error for the width 
only amounts to about 25%. This happens because of the wave function normalisation: A 
change in fs is accompanied by a corresponding change in as, both changes compensate each 
other to a large extent. Combining all uncertainties in the model parameters we estimate 
the error of the prediction to amount to about 35%. This errors applies to both the results 
presented in ( p.7|) . 

The results (p.7|) are comparable in magnitude, although in general smaller, with previous 
estimates of the perturbative contribution with two exceptions where substantially 

larger values are quoted. Simma and Wyler [Q obtain their large value by adjusting the 
perturbative contribution to the semi-leptonic decay B —>■ Deu to the data and applying 
the normalization factor obtained by that procedure also to non-leptonic B decays. Carlson 
and Milana |]6|,p!8|| use the so-called peaking approximation advocated for in In that 

approximation the i?-meson distribution amplitude is assumed to be ~ 6{x — xq) which is 
with respect to the shape of ( ^.61 ). While the peaking approximation allows one to discuss 
the qualitative features of the results in a rather simple fashion it numerically is not very 
reliable in some cases. Using it in the evaluation of the decay width within our approach, 
we find a result which is larger by a factor of 2.7 (for the AS wave) function than the value 
quoted in (|3.7|) . 

We emphasize that in the modified perturbative approach in which the transverse degrees 
of freedom and Sudakov corrections are taken into account, the soft end-point regions are 
strongly suppressed. Therefore, there is no need for an cut-off parameter in the Xi integration 
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as it is the case in the previous perturbative calculations of the B'^ decay width [Q-^]. The 
extreme sensitivity of the results to the region near xi = 1 found by these authors has 
completely disappeared in our approach. There is second advantage in our approach: we 
can make use of the standard one-loop formula for the strong coupling constant. Moreover, 
we can choose the renormalization scale as to avoid large logs from higher order perturbation 
theory (see (|2.13| )). Still the integral appearing in ( p.l9| ) is regular since the Sudakov factor 



compensates the as singularity. Actually the suppression of the end-point regions is so 
strong that the bulk of the perturbative contribution is accumulated in regions where as is 
small. In order to demonstrate that we cut off in ( p.l9|) regions where as exceeds a given 



value etc and plot Qi, as a function of ac- As can be seen from Fig. 3 50% of the result is 
already accumulated in regions where as is smaller than about 0.5. The result saturates 
for as — 0.8, there are practically no contributions from regions where as is larger than 
that value. Therefore, we consider the perturbative contribution to the decay amplitude 
as theoretically self-consistent. Note that the regions of large values of correspond to 
large quark-antiquark separations in the mesons. The use of the one-loop formula for 0:5 
(as well as the parameterization of the parton propagators) in the end-point regions may 
be questioned. However, the saturation property of the perturbative result tells us that 
as can be frozen in at a value ~ 0.8 or regularized by, say, the introduction of a gluon 
mass without changing the final results. In order to avoid the introduction of such external 
parameters we keep the standard as parameterization, being conscious of the fact that the 
actual parameterization is of no account in the soft end-point regions. 

IV. SHORT DISTANCE CORRECTIONS 

It has been shown (see, for instance, [^,0 and references therein) that leading order 
short distance corrections lead to an effective Hamiltonian at the scale /i 

7^^^ = 4%^;„ Y: a(/i)0. + 4%t;,EC.(/i)0. (4.1) 

V/ i=l.2 i=3 



replacing the weak Hamiltonian (|2.1| ) {vt = V^jV^^i = —Vu — Vc) ■ The operators Oi are given 
as follows 
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C>5 = {d^l^Lb^) J2{qpYRqp) a = (df^lf^Lb^) Y.i^c.l'^Rqp) (4.2) 

where R denotes the right-handed projection operator and q is running over the quark flavors 
being active at the scale /x. The Wilson coefficients at the scale /i = Mw are Ck{Mw) = ^fc2- 
The renormalization group evolution from the scale to /i = m^, relevant to the B decays, 



leads to the following values of the Wilson coefficients at the latter scale p6| : 
Ci = -0.229, C2 = 1.097, C3 = 0.021, 

C4 = -0.039, C5 = 0.007, Cq = -0.029. (4.3) 



At the scale /i = Mw ( [4.1|) reduces to the Hamiltonian (|2.1| ); after Fierz reordering the 
operator O2 is the one appearing in (pTl|). 

Taking into account the Hamiltonian ( |4.1| ) in the calculation of B decays into two pi- 
ous, we have to consider many contributions. There is, on the one side, the class of direct 
contributions where the d quark together with the spectator quark forms one of the pions. 
A second class contains the exchanged contributions where the d quark forms a pion to- 
gether with one of the other light quarks. Due to Fierz reordering, however, the exchanged 
contributions can be brought into the form of the contributions from the ffist class, with 
eventually different color factors and/or Wilson coefficients. Consequently, we have only to 
work out three Dirac traces and hence three integrals corresponding to the graphs shown in 
Figs. 2a, 2c and 2d. As we mentioned in Sect. 2 the graph 2b leads to a vanishing contribu- 
tion for zero pion mass. The relevant trace and the integral for graph 2a are given in ( p.7| ) 
and (|2l9| ). 



Making use of the same set of approximations as in Sect. 2, we find for the contributions 
from the graphs 2c and 2d: 

^2 = ^ J[dx]bdbb2db2K^{Xi,-h)^:^{x2,h2)'^B{x, -b) 

X exp [-S^-,{xi,b,MB, fi) - S^^{x2,b2, Mb, fi)]as{fi)[x2 + x - 1] 

xiro(V(l - a;)(l - xi) - (x - xi)x2MBb)fi^il - x)(l - xi)Mb, b, 62) (4.4) 

and 
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^3 = ^1 [dx]bdbb2db2K, i^l, {^2, h2)^B{x, -b) 

X exp [-S^-^{xi, b, Mb, //) - 5*^2(2^2, Mb, /i)] as{jj,)[xi + X2 - 2x] 



xKoiJl - xiX2 - x(2 - xi - X2)MBb)fiJ{l - x)(l - xi)Mb, b, 62). (4.5) 



In the renormalization scale expression ( p.l3| ) bi is to be replaced by 62- Obviously, the 



contributions from the penguin operators O3 and O4 lead to the same integrals, ( |2.19| ), ( [4.4| ) 
and ([4.5|) as the operators Oi and O2 respectively. The contributions from the operators O5 
and Oq, involving right-handed quarks, are either zero or lead again to the integrals ( p. 191 ), 
(0) and (^). 



There is a third class of contributions, namely those from annihilation topologies (see 
Fig. 4). In principle such topologies are generated by all six operators Oi and in each of 
these topologies a gluon can be exchanged between the quark line connecting tti and 7^2, 
and any of the other quarks, i. e., there are again four graphs contributing similar to those 
shown in Fig. 2. The contributions from the factorizing graphs are, however, zero. Only 
the contributions from the graphs where the gluon is exchanged between one of the quarks 
forming the B meson and the quark connecting tti and 112 are non-zero provided the color 
structure is appropriate. The two non-zero contributions read 

^t = ^J [dx]bdbb2db2KMM)K,{x2,h2)i'B{x, -b) 

X exp [-Sj,^{xi,b2,MB,fi) - Sj,-,{x2,b2, Mb, fi)]as{fi)[xi - 1] 

xKo(y/(1 - x)(1 - xi - X2) - X + xiX2)Mb6)/(v/(1 - xi){l - X2)Mb, b, 62) 

^l = ^j [dx]bdbb2db2K^{x^,h2)i>l^{x2,h2)i>B{x, -b) 

X exp [-S^^{xl,b2,MB,^J) - S^^{x2,b2, Mb, ^J)]as{^^)[x - xi] 



xKoi^ix - xi){x - X2)MBb)f{^{l - xi)(l - X2)Mb, b, 62). (4.6) 

Putting all together, we find the following expressions for the decay amplitudes of the 
three tttt channels 

M{B' ^ TT+TT-) = -^^v^t, {[C2 - v,/vuC43Q, 

+[Ci - vt/vuCs][n, + ^2 + ns] +[C2 - 2vt/vuC4[nt + 

M{B- 7r°7r-) = -^I^VuM%[Ci + C2] [4fi, + ^2 + ^^3)] (4.7) 
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and 

^ 7r° 7r°) = l/v^M(5° ^ TT+TT-) - MiB- n\-) (4.8) 

i. e. the amplitudes satisfy the famihar isospin relation for 5 — tttt decays [^. The contri- 
butions generated by the operators and Oq are either zero or so small (for annihilation 
topologies) that they can safely be neglected. 

Numerical evaluation of the amplitudes and the corresponding rates using the wave 
functions discussed in Sect. 3, leads to the final results for the decay widths for the three 
B decay processes shown in Table 1. The results shown in the table are evaluated with 
Vub = 0.005 and Vt = —0.0138 + zO.0050. The uncertainties due to the wave functions 
and the S-decay parameter amount to about 35%, see the discussion in Sect. 3. As an 
inspection of Table 1 reveals the QCD corrections are substantially and, in total, amount 
to a reduction of the leading term generated by the Hamiltonian (|2.1| ) of 20 — 25%. An 
exceptional role plays the decay into two uncharged pions because for that reaction 
there is no direct contribution from the Hamiltonian (|2.1| ) and the exchanged contribution 
is suppressed by color (a factor 1/3 in the amplitude). The other contributions are therefore 
relatively strong. The penguin contributions provide corrections of the order of 20% to the 
process B^ tt+tt^. Since vt is complex the CP conjugated B^ decay rates slightly differ 
from the B^ rates presented in the table (by ~ 0.2%). We stress that our aim is to estimate, 
as accurate as possible, the strength of the perturbative contributions to the B decay rates. 
For a calculation of such subtle effects as CP violations a more refined treatment of the 
penguin contributions, taking into account 0(«|) corrections, is required |3,|7|,p8|. Our 
conclusions about the strength of the perturbative contributions to the B decay rates will 
not be altered substantially by a refined treatment of the penguin contributions. Note that in 
other processes, in particular if 6 ^ s transitions are involved, the role of the penguin graphs 
may be more important (see e. g. [Q). There are also reactions in which the annihilation 
contributions are more important than in the B ^ tttt decays, e. g., in B' K^K~ where 
the direct and exchanged contributions are only generated by the penguin operators. 

Our final predictions for the branching ratios are: 

Bm'^.^.-)^(^S-\'-'' <^^' 1x10- ±35%. 

vo.oo5y y 0.44 {cz) J 
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BR{B'> 7r°7r°) = { ^ ' ^ x 10 



BR{B-^7r%-)=[^^) { -^ }xl0-'±35%. (4.9) 

We emphasize again that we favor the results obtained with the AS wave function for reasons 
explained in yet the other results obtained with the CZ wave function cannot definitely 



u.uy 




0.71 


(CZ) 


0.15 


(AS) 


0.36 


{CZ) 



be excluded for the time being. The errors quoted in ( |4.9|) do not contain the uncertainty 
in the value of Vub but it is shown explicitly in ([4.91). Since the penguin operators provide 
only small corrections also the branching ratio for the first reaction is (approximately) pro- 
portional to V^f^. Since the rate for the tt^tt" channel is much smaller than the other rates 
we refrain from quoting an error for it. 

While the decay of the B^ meson into two charmless mesons has been observed exper- 
imentally the results for individual branching ratios still suffer from large statistical and 
systematical errors. In fact only an upper limit of 2.9 x 10~^ for the process B^ tt+tt^ has 
been quoted by the CLEO collaboration. Our predictions are an order of magnitude below 
that limit, tempting us to conclude that the perturbative contributions to the B —>■ tttt 
decays are likely be too small, soft physics may still dominate these processes. Yet a definite 
conclusion cannot be drawn at present. 

For comparison we also quote the relative rates for 5° decays into two-pion states with 
isospin 2 and 

' 0.84 (AS) 

^ ' (4.10) 

0.71 (CZ) 

We observe that perturbation theory slightly favors A/ = 1/2 transitions over A/ = 3/2 
transitions in B decays. 

V. B-ir TRANSITION FORM FACTORS 

Our calculations presented in Sect. 2 and 3 implicitly contain an evaluation of the B-ir 
transition form factors at — 0- These calculations can easily be generalized to other 
values of p^. The transition form factors are defined as 
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1(7171,1 = 2 






\{7T7T,I = 




50)|2 



{n^{p^)\j]^\B'ipB)) = (pI) {pB+PiT + {vt) iPB-PiT- (5.1) 

We write the energy of the tt"*" meson in the B meson rest frame as 

E^^ = 7]Mb/2 (5.2) 

where < 77 < 1. Then the momentum transfer P2 equals M|(l — rj). Calculating the two 
form factors from the graphs 2a and 2b, we find 



iCpMl /2v^' 



vr 



2 vr V A 



n,-{i-r])n. 



1 CpMl ( 2v^' 



2 TT V A 



vt^ + {i + 7])ni 



(5.3) 



where the integral Vth is given in (|2.19| ) with the replacements: Mb by ^Mb and [2x — xi] 



by [2x — 1 + (1 — x\)ri\. The contribution from the graph 2b leads to an integral similar to 

Vti = J dxdxibdbbidbi'^l^{xi,hi)'^ b{x, -h) exp [-S^^{xi,bi, Mb, fi)] 

xa,(/i)[l - x]Koi^J{l - x)(l - xi)r]MBb)fi^il - x)r]MB, b, bi). (5.4) 



The results for the transition form factors obtained with the wave functions ( p.l| ), ( p.2| ) 
and (|3.4|) for the pion and (|3.5| ), ( p.6|) for the B meson, are shown in Fig. 5. The results 



can be trusted for values of rj between 0.2 and 1. For smaller values of r] the perturbative 
contribution is inconsistent in the sense that more than 50% of the perturbative contribution 
is built up in regions where > 0.7. The perturbative contribution to the differential decay 
rate for the semileptonic B ^ n decay does not match at r/ = 0.2 with the soft pion result 
obtained in and, integrated from rj = 0.2 to rj = 1.0 is tiny in comparison with the 
current experimental limit on the branching ratio of B^ n^l^ui. 



Li and Yu [|T^ have also calculated these form factors within the modified perturbative 
approach. While we agree with their analytical expressions for the form factors, our numer- 
ical results are about an order of magnitude smaller than theirs. We suspect that Li and Yu 
normalised the B state wrongly. 

Several other predictions for the form factor at 1] = 1 are to be found in the literature 



(BSW overlap model 0, QCD sum rules [0,^, lattice gauge theory [Q). The predicted 



values from these non-perturbative approaches range between 0.24 and 0.33. For the form 
factor at ?7 = 1 the decay width r{B^ tt+tt^) can be estimated, assuming factorization 
( p.l6| ) to hold for soft physics. One finds 



and using the quoted values for -F+(?7 = 1), the predictions for the branching ratio range 
from 0.9 x 10^^ to 1.8 x 10^^, i.e., they are rather close to the upper limit measured by the 
CLEO collaboration [Q] in contrast to the perturbative contribution. 



We have calculated the rates for B ^ tttt decays in the modified perturbative approach 
in which the transverse degrees of freedom as well as Sudakov effects, comprising gluonic 
radiative corrections, are taken into account. We believe that the perturbative contributions 
to these processes are reliably estimated: The hard scale is provided by and the Sudakov 
factor suppresses the soft end-point regions strongly so that the bulk of the perturbative 
contributions is accumulated in regions where as is sufficiently small. Therefore our esti- 
mate of the perturbative contribution can be considered as theoretically self-consistent. The 
difficulties previous authors had with the singular behaviour of the hard scattering am- 
plitude disappear in our approach, there is no need for a cut-off and correspondingly no need 
for an additional external free parameter. Also the extreme sensitivity of the perturbative 
contribution to the (soft) end-point regions disappears completely. The phenomenological 
input into our calculation, namely the mesonic wave functions - which are controlled by 
long distance physics and are, therefore, not calculable at present to a sufficient degree of 
accuracy - is fairly well constrained. The influence of various parameters and corrections 
on our results is discussed in some detail. The major uncertainties arise from the CKM 
matrix element Vub, the B-meson decay constant fs and from the pion's distribution am- 
plitude (although we favor the results obtained with the AS distribution amplitude). The 
experimental upper bound for the 5° tt+tt^ branching ratio is rather large as compared 
with the magnitude of the perturbative contribution. Although definite conclusions cannot 




(5.5) 



VI. CONCLUDING REMARKS 
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be drawn as yet, we have to be aware that B ^ tttt decays are perhaps controlled by soft 
physics. Estimates of rate asymmetries based on perturbation theory may be misleading. 

In a similar fashion as the B ^ decays we can also calculate the rates for B decays into 
other light mesons involving K, p and/or K* mesons. Depending on the particular channel 
considered the role of the operators in ( |4.1| ) other than O2, may be more important than in 
the TTTT channel. Thus, as shown in [^], the penguin graphs provide substantial corrections 
to the process B^ K^tt^. We expect similarly small perturbative contributions to the 
rates of other light meson channels. Of utmost phenomenological interest are B decays into 
channels involving D mesons since the branching ratios for these processes are typically 
two orders of magnitude larger than the upper bound for the B^ tt+tt^ branching ratio. 
This enhancement is essentially due to the fact that now the CKM matrix element He is 
involved instead of the much smaller matrix element Vub- A treatment of such channels 
within the modified perturbative approach is justified since, in contrast to the standard 
approach, the relevant parton virtualities are sufficiently large, and correspondingly is 
small enough, thanks to the transverse degrees of freedom (see (|2.6| ) and ( p.l3D ). In the light 
of our experience with the vrvr channels and backed by explorative studies we do not expect 
that the perturbative contributions to the channels involving D mesons are large enough to 
account for the experimental values. 
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APPENDIX: 

In this appendix we present a few details on the Sudakov factor used in our calculation. 
The Sudakov factor comprises those parts of gluonic radiative corrections which are not taken 
into account by the usual QCD evolution. Characteristic of it are double logs produced by 
overlapping coUinear and soft divergencies (for almost massless quarks). Examples of one- 
loop graphs responsible for radiative corrections are shown in Fig. p. In axial gauges the 
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two-particle reducible graphs, like the examples shown in Fig. ^ja, give rise to double logs, 
whereas the non-reducible graphs (Fig. Hb) only lead to single logs. The double logs from 



higher order loops can be resummed using the techniques developed by Collins et al. pQ 
The resummation results in the exponentiated one-loop corrections. 

We analyse the process under consideration, B tttt, in the B meson rest frame 
(see(|2.5|)). In this frame the pions have momenta with either a large + or a large — light-cone 



component, while the other components are zero. The same uneven kinematical situation 
holds for the quarks the pions are made up. This is exactly the kinematical situation for 
which Collins et al. [|19| derived the Sudakov factor. For exclusive processes, like B tttt, 
the Sudakov factor or more precisely the Sudakov function reads [H 



+ 



4 



q + bj 
ln(2g) + l ln(26) + l^ 



q 



,27-1 - 



-I In 

m 



2 

q — b 



\ (ln^(2g) 



-\n^{2b) 



(Al) 



where the definitions 



q = In [^Q/V2Aqcd) ; b = \n (1/6Aqcd) 



(A2) 



and 



67 10 2(3o 



9 3 27 ^ 3 ^ ' ^ 



(A3) 



are used {Pi = 102 — 38r;,//3; 7 is the Euler constant). Q/y^ is the large (+ or — ) component 
of the pion's momentum for an appropriate choice of the gauge vector. In the case of i? ^ tttt 
(see ( p.5|) ) Q equals Mb and for the B-tc transition form factor Q = rjMB- ^ is the relevant 
momentum fraction Xi or (1 — Xj), i = 1,2. Note that the coefficients in the second line 
of (|A1|) differ from previous publications [pHTl[|. In the case of the pion form factor the 
effect of these corrections to the Sudakov function results in a reduction of the perturbative 
contribution of about 2% as compared to previously reported results [0. 

The range of validity of ( |A1D for the Sudakov function is limited to not too small values of 
the transverse separation of quark and antiquark in the meson. Whenever b < ^/2/^Q (i.e.. 
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b > q), the gluonic corrections are to be considered as higher-order corrections to the hard 
scattering amphtude and, hence, are not contained in the Sudakov factor but are absorbed 
in Th- Therefore, 

s{^,b,Q) = ii h>q (A4) 

is assumed. Similarly, the complete Sudakov factor exp[— S*] is set to unity, whenever it 
exceeds unity, which is the case in the small h region. As h increases the Sudakov function 
increases as well, tending to infinity with h \/Kqcd- Consequently, the Sudakov factor 
exp[— 5] drops to zero. For hi larger than 1/ Kqcd, the true soft region, the Sudakov factor 
is zero. 

For the B meson the situation is completely different from that of the vr mesons. Due 
to the large 6-quark mass the radiative corrections only produce soft divergencies but no 
collinear ones. Consequently, double logs do not appear and, hence, the Sudakov function 
for the 6-quark is zero. For the light quark there are still both soft and collinear divergencies. 
However, because of the pronounced peak the i?-meson wave function exhibits, the bulk of 
the perturbative contribution is accumulated in regions where the momentum of the light 
quark is soft, i.e., has no large component at all. As has been shown in [|l^] the divergencies 
do not overlap in these situations. Thus, also for the light quark the Sudakov function is 
zero. 



In contrast to us, Li and Yu also associate a Sudakov function to the light quark 
contained in the B meson. Because of the shape of the i?-meson wave function only values 
of X close to Xo are relevant for which there is only a weak or even no Sudakov suppression. 
Thus, wether or not a Sudakov function for the light quark in the B meson is taken into 
account is irrelevant in praxis, the numerical values for the decay rates differ only by 5%. 
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TABLES 

TABLE L The rates for B tttt decays in units of 10^^*^ eV as predicted from the effective 



Hamiltonian (4.1). Column A: Contributions from the operator 02(^^2 = 1)) annihilations ne- 
glected. Column B: Contributions from Oi and O2, annihilations neglected. Column C: As B but 
penguin contributions included. Column D: Full result with annihilations included. 









A 


B 


C 


D 






AS 


1.04 


1.11 


0.92 


0.81 






CZ 


2.58 


2.88 


2.30 


2.21 






AS 


0.050 


0.015 


0.021 


0.034 






CZ 


0.089 


0.012 


0.043 


0.035 






AS 


0.86 


0.65 


0.65 


0.65 


B- - 




CZ 


2.03 


1.53 


1.53 


1.53 
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FIGURES 

FIG. 1. The basic graph for the decay of a 5 meson into two pions. The circle stands for the 
effective weak Hamiltonian. The quark momenta are specified. 

FIG. 2. Lowest order Feynman graphs for B decays into two pions. The internal quark and 

gluon momenta are indicated. 

FIG. 3. Saturation of the perturbative contribution: Qf, vs. ac (see text). Solid (dashed) line 
represents the perturbative contribution with (without) the Sudakov factor. 

FIG. 4. The annihilation topology. 

FIG. 5. The B-tt transition form factors vs. 77. The solid (dashed) and dash-dotted (dotted) 
lines represent the results for the form factor (F_) obtained with the AS and CZ wave functions 
respectively. 

FIG. 6. Radiative corrections. 
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